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To discuss the possibility that the superconductivities in Naa:Co02-yH20 are induced by 
electron correlation, we investigate the possible pairing symmetry based on the single-band 
Hubbard model whose dispersion of the bare energy band is obtained using the band structure 
calculation of Naa;Co02-?/H20. The superconducting transition temperature is estimated by 
solving the Eliashberg equation. In this equation, both normal and anomalous self-energies 
are calculated up to the third-order terms with respect to the Coulomb repulsion. In the case 
of spin-singlet pairing, the candidate of pairing symmetry (the maximum eigen value A^^g^ 
of Eliashberg's equation) belongs to the c?-wave(E2 representation of De group). In the case 
of spin-triplet pairing, the candidate of pairing symmetry (the maximum eigen value A^^^^ 
of Eliashberg's equation) belongs to the /y(y2_33.2)-wave (Bi representation of Dg group). It 
is found that A^^^^^ ~ '^m'ax ^''^'^ the transition temperatures of unconventional pairing state 
are estimated to be low compared with the observed temperature within our simple model. 
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Recently, Tanaka et. a/ discovered superconductivity in Na^Co02 -2/1120 {x ~ 0.35,?/ ~ 1.4) 
with the transition temperature Tc ~ 5 K.^ Na^Co02-yH20 consists of a two-dimensional tri- 
angular lattice of cobalt ions formed by a network of edge-sharing CoOe octahedra, separated 
by layers of Na ions and H2O molecules. The superconductivity in this compound has been 
intensively investigated. ^"^^ The normalized electronic specific heat data in the superconduct- 
ing state well fit the dependence. The NMR/NQR measurements have been performed 
by several groups. ^^"^^ The conclusions of the groups have been contradictory regarding the 
pairing symmetry of the superconductivity, to date. Fujimoto et al. observed an absence of 
the coherence peak just below Tc^^ and quite recently, their observation has been confirmed 
by Ishida et al}^ These results may suggest the possibility of unconventional superconductiv- 
ity in NaxCo02-yH20. In this paper, we discuss the possibility that the superconductivities 
in Naa;Co02-yH20 are induced by electron correlation. In our previous study, we investi- 
gated the superconducting instabilities of a two-dimensional triangular lattice in the Hubbard 
model with the next nearest hopping integral, by the perturbation theory with respect to 
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Coulomb repulsion U and fluctuation-exchange approximation(FLEX). We have found that 
/-wave spin-triplet pairing is the most stable in wide ranges of next nearest neighbor hop- 
ping integrals and electron number densities. The perturbation approach is sensitive to the 
dispersion of the bare energy band, by its nature. This implies that the lattice structures and 
the band filling play essential roles in the calculation of Tc, and it is important to evaluate 
Tc on the basis of the detailed electronic structure in each system. Therefore, we construct a 
two-dimensional Hubbard model whose dispersion of the bare energy band is obtained using 
the band structure calculation of Na2:Co02-yH20. Based on this Hubbard model, we estimate 
Tc by the perturbation theory with respect to Coulomb repulsion U. 

Here, we briefly mention the band structure calculation of Naa;Co02-yH20 performed 
by Nisikawa et al?^ They used the full-potential linearlized augmented plane wave (FLAPW) 
method in the local density approximation (LDA). In their calculation, the H2O molecules were 
treated as mere spacers which realistically separate two C0O2 layers in a hydrated compound 
and the Na ions were treated using the virtual crystal approximation. The obtained band 
structure has two-dimensional properties. The Fermi surfaces consist of a large "troche" -like 
hole Fermi surface around the F point and small cylindrical hole Fermi surfaces near the 
K points. The Fermi level is very close to the van Hove singularity. Co t2g orbitals take 
the main part of the density of states near the Fermi energy. The obtained Fermi surface is 
similar to the Fermi surface of nonhydrated Nao.5Co02 obtained by Singh. But the volume 
of our cylindrical hole Fermi surface near the K point is larger than the volume of Singh's 
one. There is no experimental measurements about Fermi surface of Naa;Co02-yH20 to date. 
Fermi surface is one of the important matters to consider in investigating the mechanism of 
superconductivity. So comparison of band structure calculation and experiments is desired. 

We start from quasi-particle state and use the single-band Hubbard model for discussing 
the mechanism of superconductivity. We obtain the dispersion e(k) of the bare energy band 
by using the band structure calculation of NaxCo02-yH20(a; = 0.35). e(k) can be expanded 
as e(k) = — tRexp(iR • k), where L is the set of lattice vectors of a two-dimensional 

triangular lattice. We rescale length, energy, temperature and time by a,t, j respectively 
(where a,t,kB,h are the lattice constants of the hexagonal Co plane, third hopping integral, 
Boltzmann constant and Planck constant divided by 2tt, respectively). Our Hamiltonian is 
written as 

k,(T cry^a' kik2k3k4 

where ji and U are the chemical potential and the Coulomb repulsion, respectively. The coef- 
ficients tR, are shown in Fig. 1. We calculate by solving Eliashberg's equation (Fig. 2(a)). 
In the equation, the normal self-energy and the effective interaction are obtained within the 
third-order perturbation with respect to [/(Figs. 2(b) and 2(c)). The diagrams enclosed by 
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a dashed line in Fig. 2(c) are the vertex correction terms which are not direct contributions 
from spin fluctuations. The other diagrams are included in RPA and FLEX. We call the latter 
'RPA-like diagrams' in this paper. To satisfy Luttinger's theorem, that is, the conservation 
law of particle number, we adjust the chemical potential fi using the secant method. To solve 
Eliahberg's equation by using the power method algorithm, we have to calculate the summa- 
tion over the momentum and the frequency space. Since all summations are in convolution 
forms, we can carry them out using the algorithm of the fast Fourier transformation. For 
the frequency, irrespective of the temperature, we have 1024 Matsubara frequencies. There- 
fore, we calculate throughout in the temperature region T > , where is the lower 
limit temperature for reliable numerical calculation, which is estimated to be approximately 
2.5 X lO-^{> W/{2tt X 1024) ~ 2.2 x lO'^), where W is the bandwidth; we divide a primitive 
cell into 128x128 meshes. 

The possible five candidates of unconventional pairing symmetry in our two-dimensional 
model are shown in Fig. 3. In our calculation, we obtain the maximum eigen value A^^^^ of 
spin-singlet pairing belonging to the ti-wave and the maximum eigen value spin-triplet 
pairing belonging to the fy(^y2 _^,j.2ywavc. The fy(y2-^x'^ywavc pairing state has been proposed 
in several studies of other similar models.-'^'^'^^'^^ The calculated eigenvalues Xf^g.^ ^^'^ ^raa.x 
various values of T are show in Fig. 4(a). From Fig. 4(a), we can see that when we estimate 
VF(~ ISts) ~ 0.5 eV from band structure calculation, our is estimated to be far below 1 
K, which is low compared with the observed Tc — 5K. 

To examine how the vertex corrections influence Tc, wc calculate the eigenvalues Arpa-uicc 
of Eliahberg's equation by including only RPA-like diagrams of anomalous self-energies up to 
the third order. We compare Arpa-uicc with Atopt calculated by including full diagrams of 
anomalous self-energies up to third order. The calculated eigenvalues Ap^p^ j^j^^ (for /-wave 
triplet pairing ) and AJ^pA_iike i^oi d-wave singlet pairing ) are shown in Fig. 4(b). In RPA-like 
calculation, the momentum dependence of static bare susceptibility Xo(q)O) is important for 
pairing symmetry. The calculated results of XoCQjO) are shown in Fig. 5. In this figure, we 
can observe a broad peak around M the points and a sharp peak around the F point. Prom 
Fig. 4(b), we can see that the vertex correction terms reduce the eigenvalues of spin-singlet 
pairing, on the other hand in the case of spin-triplet pairing, vertex correction terms work in 
favor of realizing the spin-triplet superconductivity. The fact that the vertex correction for a 
simple one-boson propagation is important for spin-triplet pairing, is the same as the case of 
Sr2Ru04.^^ As a result, we have A^^^^ — ^max within our simple model. 

We also calculate the self-energy S^(k,a;) using the Pade approximation (whose results 
are not presented in this paper). The w-dependence of both ReE^(k/,u;) and ImS^(k/,a;) 
(where is a Fermi momentum) near u = 0, are respectively given by ReE^(kj,(j;) oc —a; 
and ImS^(kj-,u;) oc —lo'^ This behavior is the same as that for a conventional Fermi liquid. 
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Therefore, the perturbation treatment seems to work wch in our calculation. 

In summary, we have investigated the possibility of unconventional superconductivity 
originating from the electron correlation effects, on the basis of the Hubbard model whose 
dispersion of the bare energy band is obtained using the band structure calculation of 
Naa;Co02-yH20. In the case of spin-singlet pairing, we have obtained the maximum eigen 
value A^j^x belonging to the d-wave. In the case of spin-triplet pairing, we have obtained 
the maximum eigen value Xf^^^^ belonging to the fy(^y2_^^2ywave. It is found that the vertex 
correction terms reduce eigenvalues of spin-singlet pairing, on the other hand in the case of 
spin-triplet pairing, vertex correction terms work in favor of realizing the spin-triplet super- 
conductivity. As a result, it was found that A^fj^^ — '^max the superconducting transition 
temperatures are estimated to be low compared with the observed Tc ~ 5K within our simple 
model. So we cannot decide spin-singlet pairing or spin-triplet pairing, and can not get rea- 
sonable transition temperatures of unconventional pairing state. This fact indicates that we 
have to use more a realistic model such as a multi-orbital model or the other mechanism such 
as the phonon-mediated mechanism may become dominant for the superconductivity in real 
system. 
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Fig. 1. The values of the coefficients 
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Fig. 2. (a)Eliashberg's equation. The thick Hne represents Green's function with self-energy correc- 
tion. The shaded rectangle represents the effective interaction, (b) Feynman diagrams of the normal 
self-energy up to the third order, (c) Feynman diagrams of the effective interaction up to the third 
order. Solid and dashed lines correspond to the bare Green's function and the interaction, respec- 
tively. We omit writing the diagrams given by turning the vertex correction terms in this figure 
upside down. 
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Fig. 3. The possible five candidates of unconventional pairing symmetry in our model. (a)spin- 
triplet /j^(y2_33.2)-wave(Bi representation of Dg group), (b)spin-triplet /j:(j:2_3j^2')-wave(B2 rep- 
resentation of Dg group), (c)spin-triplet p-wave(Ei representation of Dg group), (d)spin-singlet 
*a;j/(2;2_3y2)(j^2_33,2)-wave(A2 representation of Dg group), (e) spin-singlet d-wave(E2 representation 
of Dg group) 
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Fig. 4. (a)The calculated eigenvalues for various values of T (b)The calculated eigenvalues Atopt 
and ARPA-iike for various values of U 
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